The alpha-effect in a turbulent liquid-metal plane Couette flow 



G. Riidiger 

Leibniz- Institut fur Astrophysik Potsdam, An der Sternwarte 16, D- 14482 Potsdam, Germany and 
Helmholtz-Zentrum Dresden-Rossendorf, P.O. Box 510119, D-01314 Dresden, Germanflj 



o 

(N 
C 

m 



in 
6 

in 

o3 



> 

in 
q 

o 



X 



A. Brandenburg 

Nordita, AlbaNova University Center, Roslagstullsbacken 23, SE- 10691 Stockholm, Sweden and 
Department of Astronomy, Stockholm University, SE- 10691 Stockholm, Sweden 

(Dated: January 4, 2012) 

The electromotive force in a plane Couette flow of conducting fluids under the influence of a large- 
scale magnetic field is calculated for driven turbulence. In case of a vertically stratified intensity 
of the turbulence the resulting a-tensor in presence of a horizontal shear also possesses diagonal 
elements known as the so-called a-effect. The possibility to measure these components is probed 
by both quasilinear theory and nonlinear numerical simulations. For all magnetic Prandtl numbers 
< 1 we find the a-component along the flow in the high-conductivity limit exceeding the other 
components. In this limit the cv's run linear with the magnetic Reynolds number Rm while in the 
low-conductivity limit they run with the product Rm ■ Re. For the small magnetic Prandtl numbers 
of liquid sodium a value for the electromotive force across the channel of order 0.05. ..0.5 mV for 
magnetic fields of about 1 kGauss and velocity fluctuations of about 1 m/s is found for a channel 
of (say) 50 cm height independent of its width. Such values should be appropriate to be realized in 
a future laboratory experiment. 

PACS numbers: 47.27.Pa, 47.65. +a, 95.30.Qd 



I. INTRODUCTION 

The electrodynamics in turbulent conducting fluids 
provides the standard explanation of the existence of 
magnetic fields of the cosmic objects. The excitation of 
fields results from the interplay of two different elemen- 
tary processes. It is known that the turbulent motions 
are reducing large-scale electrical currents by inducing 
an electromotive force (EMF) opposite to the direction 
of the current. One can write 



-HoVtJ, 



(1) 



where rfx is the turbulent magnetic diffusivity. In stel- 
lar convection zones its value exceeds the molecular (or 
'microscopic') value of the magnetic resistivity by many 
orders of magnitudes. 

In this paper the EMF is derived for a turbulent fluid in 
the presence of a shear flow U and a uniform background 
field B. The fluctuating flow components are denoted by 
u and the fluctuating field components are denoted by b. 

For a turbulent dynamo the enhanced dissipation must 
be overcome by an induction process which does not run 
with the electric current. One also knows that under the 
influence of global rotation and a uniform magnetic field, 
anisotropic turbulence produces an EMF parallel to the 
field, i.e. 



£ = aB 



(2) 



Here, a is a pseudo-scalar formed by the rotation vector 
17 and the anisotropy direction g with a oc g ■ f2. Often 
as the anisotropy direction is the gradient of the density 
stratification in the fluid is used, but it is also possible 
that an intensity gradient close to rigid boundaries forms 
the preferred direction. 

The resulting dynamo equation in rotating and strati- 
fied plasma is 



dB 



dt 



— = rot [aB-(r) + r] T )mtB ) . 



(3) 
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which indeed has non-decaying solutions if the a-effect 
exceeds a critical value ('a 2 -dynamo'). 

In many of papers the presented concept of a turbu- 
lent dynamo has been applied to planets, stars, accretion 
disks, galaxies and galaxy clusters (many references given 
by Only very few papers, however, deal with an ex- 
perimental confirmation of the validity of relations ([T]) 
and ([2]) in the laboratory. This is surprising insofar as 
Eq. © as a simple consequence of Eqs. and seems 
to summarize all our knowledge about the electrodynam- 
ics in rotating turbulent fluid conductors. Generally, the 
validity of Eq. (P) is not seriously doubted. However, the 
existing laboratory experiments report an increase of the 
molecular value of r\ by only a few percent. This is a 
surprising situation as the eddy-concept of the effective 
dissipation in turbulent media governs much of the cos- 
mical physics from the climate research, geophysics up to 
the theory of star formation and quasars. 

An even more dramatic situation holds with respect to 
the a-effect. There are one or two existing experiments 
on the basis of the idea that the a-effect is essentially 
a measure of the swirl of the fluid. It has indeed been 
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shown that a fluid with imposed helicity (imposed by 
rigid, swirling channels), produces an EMF in the direc- 
tion of an imposed field; see 0-0] • It is not yet shown, 
however, that a rotating fluid with helicity which is not 
imposed (but results from the global rotation of the fluid) 
leads to an observable a-effect. In nature, helicity is usu- 
ally due to the interaction of rotating turbulence with 
density stratification. The present proposal also aims in 
the direction to probe the possibilities of such a more 
rigorous a-effect experiment. The many difficulties for 
such an experiment which we shall demonstrate make it 
understandable that an experiment without prescribed 
helicity not yet exists in the theory and in the MHD lab- 
oratory. 

It is easy to formulate the overall problem of all a- 
effect experiments. With Eq. ([2]) the potential difference 
between the endplates of the container in direction of the 
mean magnetic fields is 



A<P = UdynBH 



(4) 



with H as the distance between the endplates. Hence 
for (say) H ~ 100 cm and (say) B ~ 1000 Gauss the 
potential difference is A4> — ad yn | in mV and cm/s, re- 
spectively. The maximum a-value is of the order it rms 



hence < w rms in mV. For u T 



1 cm/s the induced 



potential difference, therefore, is maximally 1 mV. A con- 
tainer of 5 cm radius rotating with 1 Hz has a linear 
outer velocity of more than 30 cm/s so that it rms ~ 1 
cm/s might be considered as a conservative estimate. We 
find as a necessary condition for any a-experiment that 
one must be able to measure potential differences smaller 
than a few mV. The a-experiment in Riga worked with 
B ~ 1 kGauss and velocities of order m/s so that the A<P 
exceeded 10 mV. This experiment, however, used a pre- 
scribed helical geometry to mimic the symmetry breaking 
between left and right helicity @- 

If the basic rotation is not uniform then its shear 
induces toroidal magnetic fields so that for sufficiently 
strong shear even the a-effect can be rather small ('af2- 
dynamo'). Also turbulence subject to a plane shear flow 
(without rotation!) of a liquid metal is able to work as 
a dynamo if the turbulence intensity is stratified in the 
direction(orthogonal to the shear flow plane, see [5J.A ba- 
sic rotation may thus not be the only flow system whose 
influence enables the turbulence to generate global fields. 

In the present paper a plane Couette flow is considered 
to analyze the characteristic issues of the corresponding 
a effect and to design a possible experiment to measure 
its amplitude. 



II. THE MEAN ELECTROMOTIVE FORCE 

Consider a plane shear flow with uniform vorticity in 
the vertical z direction, i.e. 



The shear flow may exist in a turbulence field which docs 
not possess any other anisotropy than that induced by 
the shear ([5]) itself. The one-point correlation tensor is 



Ui{x, t)Uj(x, t). 



(0) 



The correlation tensor may be constructed by a perturba- 
tion method. The fluctuating velocity field is represented 
by a series expansion, 



,(°) 



(7) 



where the upper index shows the order of the contribu- 
tions in terms of the mean shear flow. 

The zero-order term represents the 'original' isotropic 
turbulence which is assumed as not yet influenced by the 
shear. The spectral tensor for the original turbulence is 



Q 



(o) _ E(k,u) 
Wirk 2 



Si 



(8) 



where the positive-definite spectrum E gives the intensity 
of isotropic fluctuations, i.e. 



/.(or = 



OOOG 

J j E(k 7 uj) dk duj. 
o o 



(9) 



For analytical calculations the one-parametric spectrum 

(10) 



2 w 

E{k,u>) = ——E(k), 



Uy = SX. 



(5) 



7T U! 2 + W 2 

can be used which provides E oc S(lu) in the limit of w — ¥ 
and leads to white noise for large w. The correlation 
time of the turbulence is defined by r corr = 1/w. 

By definition, the magnetic-diffusivity tensor relates 
the mean electromotive force (TT|) to the gradients of the 
mean magnetic field via the relation Si = r\ijkBj,k- This 
tensor for an isotropic turbulence influenced by a mean 
field shear flow §5§ has been constructed up to the first 
order in the shear @. It has been shown that the 
combination of the shear and the shear-induced parts of 
the magnetic-diffusion tensor are not able to operate as 
a dynamo. 

It has also been shown that the shear in combination 
with a stratified turbulence provides a helical turbulence 
which leads in Eq. ^ to an a-effect [5| . Here a must be a 
pseudotensor so that an e-tensor has to appear in the a- 
coefficients. The construction of the EMF £~i — tijkUjbk 
is the only possibility for the e-tensor to appear. There- 
fore, the subscript of Si is always a subscript of the e- 
tensor. As the e-tensor is of 3 rd rank an inhomogeneity of 
turbulence with the stratification vector, g = V log u 2 ms 

with u rms = Vu^, must also be present for the a-effect 
to exist. If the shear flow is included to its first order, 
the general structure of the a-tensor is 

atij = -ye ijk gk + {ot\£ikiUj,k + a 2 eikiUk,])gi + 

+ a 3 e lk igjUi,k + aA^kjUi,k9i + a 5 e ijk Uk,igi(a) 
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If the stratification is along the vertical z-axis it follows 
from |QT]) that 

Olxx = OL2QzS = a x S, 

a yy — -a\g z S = a y S, 

at zz = a 3 g z S = a z S 

a xy = -a V x = ig z = r. (12) 

The turbulent pumping is described by a xy . The 
anisotropy of the a-tensor is described by the difference 
between a x and a y . In the used geometry the a yy (i.e. 
the coefficient ot\ defined below) plays the main role in all 
cosmical applications while in the proposed experiment in 
a turbulent shear flow the coefficient ui is probed which 
produces the EMF perpendicular to the flow. 
The coefficients of (fTTj) read 



7= 6 



r)k 2 E(k,uj) 
uj 2 + rj 2 k 4 



dfc da; 



(13) 



for the pumping term and 

oooo 

a n = J J A n E(k,Lj)dk du, u=l,2,...5, (14) 
o o 

for the a-effect with 

Avrfk® + 2u?r\ (y + if) k 4 



A 2 



A, 



15 (w 2 + v 2 k 4 ) (uj 2 +r] 2 k 4 y 
rfk 4 (r] 2 k 4 - 3uj 2 ) 
15 {uj 2 + r] 2 k 4 ) 3 ' 

rfv 3, (4?7 — hv) k 12 
~60 (w 2 + v 2 k 4 f (lu 2 + T) 2 k 4 f 
u?v (28?7 3 - 4r; 2 ^ + Ylr)v 2 + hv z ) k 8 
60 (uj 2 + v 2 k 4 f {u 2 + 7] 2 k 4 ) 2 
(jj 4 rj (if + 36v) k 4 — 5w 6 
60 (w 2 + v 2 k 4 f (uj 2 + t] 2 k 4 ) 2 
8uj 2 k 4 (u} 2 + vr\k 4 )ri(v + if) 
~ 30(uj 2 + v 2 k 4 ) 2 (uj 2 + r] 2 k 4 ) 2 ~ 
3?7 4 fc 8 -24w 2 r7 2 fc 4 + 5cj 4 
30(w 2 + rj 2 k 4 ) 3 



(15) 



for the kernels. Only the terms occurring in (1121) have 
been given. Note that A3 for v — > exists 



A, = 



11t7 4 A: 8 - 16wVfc 4 



5w 4 



30(o; 2 + ry 2 fc 4 ) 3 
For small Pm the expression for A\ simplifies to 

1 1 \ Avrfk 9, 



(16) 



A : 



15 (uj 2 +r] 2 k 4 f \u: 2 + v 2 kA 
ifk 4 (3?? 2 fc 4 -w 2 )l 

UJ 2 + n 2 k 4 | 



(17) 



The first expression on the RHS for v — > forms a (5- 
f unction, so that 



2tt 
15ui 



(i 



E(k,0) 
rjk 2 



dfc 



with 



Rm 



V 



3?7t 
V 



(18) 



(19) 



The missing terms in (|18[) result to be of the same order. 

The white-noise approximation mimics the high- 
conductivity limit which holds for cosmical applications. 
In this approach the spectrum does not depend on the 
frequency u up to a maximum value w max from which one 
the power spectrum vanishes. It is a turbulence model 
with very short correlation time, i.e. r corr ~ l/oj max . One 
finds from (|17p after integration 



n fE(k,0) 



677 



fc 2 



dfc. 



(20) 



The resulting integral can be written as 7" C orr^corr u rms 
that 



so 



ai = — Rm 
6 



I 2 

w corr 1 



(21) 



which hardly differs from the result (fT5)l . As it must the 
factor 7r/6 also appears in Fig. Q] (middle, bottom) as the 
value of Ji for t = 0. For Pm = 1 the corresponding 
expression is 



ol\ = — Rm 
15 



I 2 

COIT ' 



(22) 



Note that the value for small Pm exceeds the value for 
Pm = 1. For given r\ smaller v lead to higher values of 
the EMF. 

The full expression might be written in the form 



a\ ~ I\ Rm 1, 



2 

corr J 



(23) 



where I\ is numerically given in Fig.[T] It is the result of a 
numerical integration under use of a spectral function of 
the form exp(— T 2 a> 2 ) with the dimensionless parameter 



t = i]k 2 T c . 



1/Rm. 



(24) 



While small r represent the high-conductivity limit the 
larger r represent the low-conductivity limit. Obviously, 
the high-conductivity limit can be modeled with a white- 
noise spectrum (Eq. (|10[) with w — > 00) while the low- 
conductivity limit can be treated by use of the 5-function 
(Eq. flU with w-tO). It follows from (fl"5|h 



1 / 4 

a>i = — H 

15 V Pm 



Rm 2 it 



(25) 
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FIG. 1: The numerical values of the coefficients Ji, I2 and 
J 3 . Top: Pm = 1, middle: Pm = ICT 5 , bottom: Pm = ICT 6 . 
Note that in the high-conductivity limit (r — > 0, white noise 
spectrum) the coefficients I\ and I2 for Pm <C 1 exceed the 
values for Pm = 1. Laboratory experiments with turbulent 
flows of liquid metals work with r ~ 1...10. 

so that Pm = 1 leads to a± = |Rm 2 £^ OI1 . and for small 
magnetic Prandtl number one finds ot\ — j^Re Rm ^ orr , 
with Re = Rm/Pm as the Reynolds number. Hence, in 
the low-conductivity (small Rm) limit the a effect runs 
with Re • Rm while in the high-conductivity limit (large 
Rm) the a effect runs with Rm. Very similar expressions 
also occur if the shear flow is formally replaced by a basic 



FIG. 2: Ii for various Pm. Note that for given magnetic dif- 
fusivity rj fluids with much smaller viscosity (like sodium and 
gallium) are even better qualified for laboratory experiments. 

rotation (Riidiger 1978). 

Figure [5] demonstrates the behavior of the numerical 
integrals for decreasing Pm and growing r. The integrals 
are running with (1 + Pm r) _1 so that for Pm t> 1 the 
relation a\ oc Re Rm indeed results. Again, for given 
magnetic diffusivity smaller values of the viscosity en- 
hance (!) the resulting cv-effect. 

In experiments with liquid metals the Rm is of order 
0.1. ..1 so that t ~ 1...10. In this domain and for small 
Pm the coefficient 1\ hardly changes with t (see Fig. [I}. 
Its approximate value is 0.5 which is very close to the 
value 7r/6 valid in the high-conductivity limit. The rea- 
son is that for small Pm the transition of 1\ to the low- 
conductivity limit only happens at rather high values of 

T(Fi g .rj. 

Similar calculations for a 2 lead to a 2 — —I 2 Rm ^ orr , 
with I 2 also plotted in Fig. [T] For medium r and small 
Pm we find 1% — I\, i.e. 

a 2 ~ -0.5 Rm £* „. (26) 

In the low-conductivity limit we have 

a 2 = -i(A_ 5 ) Rm2 , Lri (2?) 

which changes its sign at Pm = 0.8. For small Pm it can 
also be written as 

a 2 = -^RmRe £ 2 cm ., (28) 

In summary, the plots in Fig. Q] reveal an important in- 
fluence of the viscosity value for given diffusivity value. 
For small Pm in the low-conductivity limit the ratio of I 2 
to 1\ approaches unity while it is very small for Pm = 1. 
The magnetic Prandtl number is highly different for liq- 
uid metal and for numerical simulations. Figure [T] shows 
that in numerical simulations the a xx should be much 
smaller than the a vv what is not true for laboratory con- 
ditions with their small Pm. 
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The «3 yields «3 ~ — I3 Rm £ 2 orr with J 3 of about 0.1 
for small Pm (see Fig. [T]). The horizontal components 
of the a-tensor are thus of the same order of magnitude 
but the vertical component appears to be smaller. Note 
that for Pm = 1 the a\ strongly exceeds the cti but it 
does not for small Pm. The signs of all components are 
identical. 

The pumping term does not depend on the shear. After 
(1131) for small 77 it does not run with I/77, so that it is 
simply 
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2 >T C , 



77 Rm. 



(29) 



It can only be measured if the external field Bq lies in 
the shear plane. 



III. NUMERICAL SIMULATIONS 

It is straightforward to verify the existence of an 
a-effect using numerical simulations of non-uniformly 
forced turbulent shear flows in Cartesian coordinates. 
We perform simulations in a cubic domain of size 
so the minimal wavenumber is k = k± = 2tt / ' L. We solve 
the governing equations of compressible hydrodynamics 
with an isothermal equation of state with constant sound 
speed c s . The set of governing equations is 

^- = SU x y-c 2 s V\np+f + p- 1 V-2pvS, (30) 



Dt 



Dlnp 
D/j 



-V-U, 



(31) 



where D/D/j = d / dt + U ■ V + Sx d/dy is the advective 
derivative with respect to the full velocity field (includ- 
ing the shear flow), U is the departure from the mean 
shear flow, and = \{diUj + djUi) — V • U is the 
trace-less rate of strain matrix. The flow is driven by 
a random forcing function consisting of nonhelical waves 
with wavenumbers whose modulus lies in a narrow band 
around an average wave number kf = 5ki 0]. We arrange 
the amplitude of the forcing function such that the rms 
velocity increases with height while the maximum Mach 
number remains below 0.1, so the effects of compressibil- 
ity are negligible. 

We use the kinematic test-field method [9J in the Carte- 
sian implementation fioj ] to compute from the simula- 
tions simultaneously the relevant components of the a 
effect and turbulent diffusivity tensors, a y - and rjij, We 
do this by solving an additional set of equations govern- 
ing the departure of the magnetic field from a set of given 
mean fields. This mean field is referred to as the test field 
and is marked by the superscript T. For each test field 

B , we find a corresponding fluctuation 6 T = V x a T by 
solving the inhomogeneous equation for the correspond- 
ing vector potential a T , 



Pa 1 
Vt 



-Sa^x + Uxb 1 



MX B 



u x b 1 



-r/V 2 a T , 
(32) 



where T>/T>t = d/dt + Sx d/dy is the advective deriva- 
tive with respect to the imposed shear flow only, and 
(« x 6 T ) = u x b T — u x b T is the fluctuating part of 
u x b T . We compute the corresponding mean electromo- 



tive force, £ = u x b T , which is then related to B and 

— t — T 

its curl, [XqJ =VxB , via 



_t — T 



1 T 



(33) 



We use 4 different test fields with x or y components 
being proportional to sin kz or cosfcz. The x and y com- 
ponents of Eq. ([3"3"]l constitute then 8 equations for the 4 
relevant components of a,ij{z,t) and r]ij{z,t). 

We adopt periodic boundary conditions in the y direc- 
tion, shearing-periodic boundary conditions in the x di- 
rection, and stress-free perfect conductor boundary con- 
ditions in the z direction, i.e., 



d z 



d z a, 



0. 



(34) 



A numerical resolution of 64 3 mesh points was found to 
be sufficient. The Pencil Code 18j has been used for 
all calculations. 

Simulations are performed for a number of different 
parameter combinations. The quantities S and g z are 
positive in the calculations presented here, i.e. the ba- 
sic velocity ([5]) grows in the positive x direction while 
the turbulence intensity grows in the positive z direc- 
tion. To make contact with laboratory experiments, we 
focus here on the case of low conductivity and choose 
Rm = M rms /i)fcf = 0.2, which is consistent with our defi- 
nition of Eq. (|19p with a Strouhal number of unity, i.e., 



' corr ^rms 



kf = 1. 



As in earlier work using fully helical turbulence, we 
present the components of a y and 77.^ in normalized form 
in terms of ao = u lms /3 and 77x0 = Urms/3fcf, respec- 
tively. The shear of the background flow is normalized 
with the speed of sound, i.e. 



S = SC ac fci 



2nsu r 



M&D 



(35) 



with the Mach number Ma = ii, ms /c ac . Here, D (= L) 
is the distance between the walls, i.e. the width of the 
domain in the x direction. In the simulations we work 
with s = 0.2 and Ma = 0.05. Hence, 



yy 
ao 



-12 7i s - 



Rm 



D 2 Ma 



(36) 



l yy\ 



Following Eqs. p^|) . both the streamwise and the span- 
wise a-tensor components, a yy and a xx , should be neg- 
ative. As the simulations are done for Pm = 1 the \a v 
should strongly exceed the value of la^l. 

Figure shows these predictions as fulfilled. They 
hold for Rm = 0.2 and kf/ki = 5 so that 10 cells exist in 
the vertical direction. From Eq. (|3"o| one finds a vy /a ~ 
—0.01. Both the a-tensor components are negative and 
\a yy \ dominates. Its amplitude is indeed 0.01. For Pm = 
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FIG. 3: Simulations with positive shear (s = 0.2): The numerical values for a (top) and the two shear-current terms r] xy and 
r/yx (bottom). 



1 the I-x is strongly reduced relative to I\ so that the small 
amplitude of a xx in Fig. [3^, becomes understandable. For 
small magnetic Prandtl number this reduction does not 
exist (see Fig. [1]). The simulations also confirm that the 
signs of the diagonal components of the a tensor only 
depend on the sign of the product g z S, as formulated in 
the relations (fl"2"j) . 

Next, we look at the off-diagonal components of o^; 
see Fig. Note that a yx « — a xy , which corresponds 
to a turbulent pumping velocity in the z direction. It is 
negative for k\z < 2.5, corresponding to downward trans- 
port, i.e., down the gradient of the turbulent intensity, as 
expected 0j- Near the top of the domain, the gradient 
of the turbulent intensity is reversed and so is the sign of 
the pumping velocity, ctyx- 



IV. SHEAR— CURRENT EFFECTS 

Let us now discuss the components of the rjij tensor. 
Before we turn to the shear-current effect, we note that 
the two relevant diagonal components are nearly equal, 
i.e., -q xx « i] yy (Fig. |3t). This high degree if isotropy of 
turbulent diffusion in the xy plane has already been in 



earlier measurements of turbulent diffusivity in turbulent 
shear flows [TollTTj. The maximum of i]t = {i]xx+Vyy)/2 
is about 20% of the reference value t^to, which agrees 
with the fact that for Rm <§; 1, ??t/?7to ~ R- m [HI- 

Let us turn to the off-diagonal components. For non- 
stratified shear flows one finds 



7]xy — S, Tjy X — Tjy S 7 



(37) 



so that the 77-terms depend on the sign and the amplitude 
of S. In a quasilinear theory (SOCA) the X] y proves to 
be negative-definite [5j . For positive shear the coefficient 
i] yxi therefore, is expected to be negative. This result has 
also been confirmed for Rm < 200 [1(| ■ The possibility of 
dynamo action in these nonhelical shear flows [T3, Il3| is 
therefore not explained by the shear-current effect 11411 . 
but instead by the incoherent alpha-shear dynamo [15| . 

The same dependence of the sign of r\ yx on that of S 
holds also for rj xyj but only for Pm of order unity and 
in the low conductivity limit (small Rm). Both the con- 
ditions arc fulfilled for our simulations. Experiments of 
liquid metals, however, concern small magnetic Prandtl 
numbers for which the rj x results as positive. As simu- 
lations for such small Pm are not positive, only experi- 
ments can finally provide the sign of r\ x . The numerical 
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simulations with positive shear and Pm = 1 indeed pro- 
duce the expected negative values for both r\ xy and r\ yx 

(Fig. EH). 



SHEAR-FLOW ELECTRODYNAMICS 



same. Note that the statement bases on the assumption 
^cott/H ~ const, in both cases. Note also that turbulence 
which is stratified in the z direction, there is no pumping 
if also the magnetic field points in the spanwise direction. 



The electromotive force across a channel after (fl2t is 



£x = a 2 g z SB 



(38) 



so that the potential difference S'Px between the walls 
with their distance D is 

5$ x = a 2 g z SDB ~ -0.5Bjnf cmi g z UB 

~ -Rm A 2 HUBq (39) 

with g z = —2/H as the vertical scale height of the tur- 
bulence stratification and the ratio A = £ con -/H. The 
amplitude of the mean shear flow is U. Note that the 
width D of the channel does not appear in (|39|) . Hence, 



8$ x [mV] ~ 10 RmA 2 



H 



10 cm 



U 
m/s 



Bn 



kGauss 



(40) 



so that with (say) A ~ 0.1 and a channel of 50 cm height, 
a shear flow of 1 m/s subject to a magnetic field of 1000 
Gauss leads to a potential difference of 



5<P X ~ 0.5 Rm [mV]. 



(41) 



For the (maximal) value of Rm ~ 0.1.. .1 (u rms ~ lm/s 
and ^corr — 5 cm) the measurements may provide a po- 
tential difference of 0.05. ..0.5 mV. The numbers are quite 
similar to those of the Riga experiment 0, 0] • The basic 
difference is that in a shear flow experiment, the helicity 
is not prescribed but it is naturally produced by the inter- 
action of the (stratified) turbulence with the background 
shear. 

If the external field is perpendicular to the shear plane, 
then the result for the potential difference in the z direc- 
tion results as 



5<P Z ~ 5<P X 



Din 



(42) 



After Fig. [T]the ratio I3/I2 is of the order of 0.1 so that 
for D ~ H the spanwise potential difference is smaller 
than the same value in the perpendicular x direction. 
For narrow channels with D ~ 0.1H both values are the 



VI. CONCLUSIONS 

Laboratory studies of homogeneous dynamos are still 
in their infancy. The only working dynamo where the 
flow pattern is not stronhgly constrained by pipes or con- 
tainer walls is the experiment in Cadarache [16| where, 
however, the effects of soft iron play an important and 
not well understood role. No contact with the a-effect 
in mean-field theory has been possible in that case. The 
present proposal of measuring the a-effect in an uncon- 
strained turbulent flow would therefore be a major step 
forward. In such an experiment, the pseudo-scalar neces- 
sary for producing helicity comes from the stratification 
of turbulent intensity giving rise to a polar vector and the 
vorticity associated with the shear flow giving rise to an 
axial vector. It has been shown, therefore, that the basic 
effects of the theory of the turbulent dynamo which are 
usually considered as special properties of rotating fluids 
can also be found for the plane-shear flow. 

The present work yields a detailed prediction about 
the sign and magnitude of the tensor components of onj 
and rjij . It may motivate the construction of a suitable 
experiment using liquid metals to achieve a measurable 
a-effect. The needed vertical stratification of the turbu- 
lence intensity must experimentally imitated using grids 
with nonuniform mesh sizes and/or by walls of increas- 
ing/decreasing roughness of the walls in vertical direc- 
tion. 

We have shown that the values of a are small compared 
with naive estimates. In fact, it may not be possible that 
such flows could produce a supercritical dynamo in the 
conceivable future. Nevertheless, even in the subcritical 
case, an a-effect should be measurable, which would thus 
open the possibility of detailed comparisons between the- 
ory, simulations, and experiments. Once such a compar- 
ison is possible, there will be more details that should be 
investigated. One of them concerns the modifications of 
the results in the presence of imperfect scale separation 
in space and time. For oscillatory dynamos, this effect 
can significantly lower the excitation conditions for the 
dynamo compared to standard mean-field estimates [I?) ■ 
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